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ABSTRACT 

Recent X-ray variability studies on a large sample of Active Galactic Nuclei (AGN) suggest that the log of 
the square of the fractional rms variability amplitude, <r r 2 ms , seems to correlate with the log of the black-hole 
mass, Mbh> with larger black holes being less variable for a fixed time interval. The fact that the rms amplitude 
can be easily obtained from the observed X-ray light curves has motivated the theoretical modeling of the a^ ms - 
Mbh correlation with the aim of constraining AGN masses based on X-ray variability. A viable approach to 
addressing this problem is to assume an underlying power spectral density with a suitable mass dependence, 
derive the functional form of the of ms - Mbh correlation for a given sampling pattern, and investigate whether 
the result is consistent with the observations. Inspired by the similarities shared by the timing properties 
of AGN and X-ray binaries, previous studies have explored model power spectral densities characterized by 
broken power laws. For simplicity these studies have in general ignored the distorting effects that the particular 
sampling pattern used to obtain the observations imprints in the observed power spectral density. With the 
advent of RXTE, however, it has been shown that the band-limited noise spectra of X-ray binaries are not 
broken power laws but can often be described in terms of a small set of broad Lorentzians, with different 
amplitudes and centroid frequencies. Motivated by the latest timing results from X-ray binaries, we propose 
that AGN broad-band noise spectra consist of a small number of Lorentzian components. This assumption 
allows, for the first time, to fully account for sampling effects in theoretical models of X-ray variability in an 
analytic manner. We show that, neglecting sampling effects when deriving the fractional rms from the model 
power spectral density can lead to underestimating it by a factor of up to 80% with respect to its true value for 
the typical sampling patterns used to monitor AGN. We discuss the implications of our results for the derivation 
of AGN masses using theoretical models of the of ms - Mbh correlation. 

Subject headings: accretion disks — black hole physics — galaxies:active, nuclei — X-rays:binaries, galaxies 



1. INTRODUCTION 

Most of the known accreting compact objects emitting in 
the X-ray band, both in our and other galaxies, show large 
variations in their fluxes (up to a factor of a few or even more 
in some cases) over several decades in frequency. X-ray bi- 
naries, hosting neutron stars and black holes with masses be- 
tween 1M and 10 M , are known to be variable on time 
scales ranging fro m milliseconds up to several days (see 
Ivan der Klisl 1200 5. for a recent review). The X-ray emis- 
sion from Active Galactic Nuclei (AGN), on the other end 
of the mass spectrum, with masses in the range 10 6 M Q to 
10 9 Mq, fluctuates on tim e scales ranging from a few m inutes 
up to a few years ("e.g.. iMa rshail., Warwick. & Poundsl [l981t 
Mushotzk v. Done. & Pou nds 1993). The origin of this strong 
variability is currently not understood. 

From the theoretical point of view, the main goal of vari- 
ability studies is to understand the physical processes that 
modulate the emission of high-energy radiation and, in par- 
ticular, how the fundamental properties of the central com- 
pact object (mass, spin, surface or lack thereof) determine the 
spectrum of observed frequencies. From the phenomenolog- 
ical point of view, observations provide a powerful tool for 
constraining mathematical models which, in turn, provide a 
guide in the investigation of physical models that can account 
for the observed variability. The studies carried out over the 
last decade with the Advanced Satellite for Cosmology and 
Astrophysics (ASCA) and the Rossi X-ray Timing Explorer 
(RXTE) and, more recently, with the Chandra X-ray Obser- 



vatory and XMM-Newton have revolutionized the timing phe- 
nomenology of both X-ray binaries and AGN. 

In the case of galactic sources, because of the typical 
timescales involved and thanks to the unprecedented tim- 
ing capabilities of RXTE, it has become possible to derive 
the power spectral density (hereafter the power spectrum) 
of a larg e number of a ccreting binaries with exquisite detail 
(see e.g., [van der Kli a 120051) . This made possible an impor- 
tant step forward toward understanding the similarities and 
differences exhib ited in the power spec t ra of neutron stars 
and black holes (Mivamoto et al. 1994; van derKlis 1994; 
Sunvaev & Revnivtsev 2000; Belloni. Psaltis. & van der Klis 
120021) 

Recently, it has become evident that the most salient 
timing features present in many X-ray binaries can be 
well described by a small number of Lorentzians with 
different centroid f requencies, widt hs, and amp litudes 
JNowakll2000l iBelloni. Psaltis. & van der Klisl 12002). This 
phenomenological framework has provided a solid ground 
for studying a number of tight correlations exhibited by 
the parameters characterizing the different Lorentzian 
components. Surprisingly, these correlations are not only 
present in any given source but are also maintained across 
sources (both neutron stars and b lack holes) over several 
decades in frequency space ( Wiinands & van der Klis 
[i999t iPsaltis. Belloni & van der Klisl 119991: 

Belloni, Psaltis, & van der Klis 2002). Some of these 
correlations might even be present in white-dwarf systems, 
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extending the range of their validi ty to lower frequ e ncies by 
two o rders of magnitude ( Warner & Woudtj 120021 [Mauche 
2002). The study, and eventual understanding, of these 
correlations offers one of the most promising avenues 
for constraining theoretical models of X-ray variability in 
accreting binaries. 

The lower observed fluxes and longer timescales associated 
with AGN have not made possible a comparable progress in 
our understanding of their timing properties. Both of these as- 
pects affect the characterization of the power spectrum from 
the data at high and low frequencies, but also, through sam- 
pling effects, within the accessible frequency range. If there is 
significant variability below or above the lowest and/or high- 
est frequencies probed by the observations, then the estima- 
tion of the power spectrum from the data suffers fr om the ef- 
fects of red noise le akage and aliasing (see, e.g.. i DeemingI 
[T975llPriestlevl1 989tly^nder Kli Jl 989tlPress et al J1 99%L 

In the last few years, however, it has become clear that AGN 
power spectra do present characteri stic frequencies ("breaks") 
on the scale of hours to months (Edelson & Nandral ll999i 
IUttlevetalJl2002t iMarkowitz 'et al] 120031) . These frequen- 
cies have been linked to similar characteristic frequencies de- 
tected in early observatio ns of a number of galacti c black- 
hole candidates (see, e.g., Belloni & Hasinger 1990). More- 
over, there is a growing body of evidence suggesting that the 
similarities shared by the X-ray variability of AGN and X- 
ray binaries extends beyond the overall shape of their power 
spectra. These similarities encompass a strong linear correla- 
tion between the rms variability amplitude and the X-ray flux, 
time scale-dependent lags between the hard and soft X-ray 
bands, and a similar energy-dependence of the power spec- 
trum shape with hig her energies showing flatt er slopes above 
the frequency break ( Uttlev & McHardv 2004, and references 
therein). The remarkable similarities shared by the global tim- 
ing properties of objects with masses that differ so widely and 
the fact that the characteristic dynamical frequencies observed 
in AGN seem to correlate with mass, has encouraged the use 
of X-ray variability as a mean for estimating supermassive 
black-holes masses. Two different approaches are in use. 

For the handful of AGN for which there are good available 
data covering a long period of time, it is possible to iden- 
tify one (and in a few cases two) characteristic frequency(ies) 
in their power spectrum. Assuming that the same charac- 
teristic frequency can be identified in a galactic black hole 
of known mass (usually Cyg X-l) and that the characteristic 
break frequency scales in versely with mass, the mass of sev- 
eral AGN can be inferred jlJttlev et alJ2002tlMarkowitz et al.1 
l2003llMcHardv et al 120041120051) . 

As it is most often the case, however, the sparsity of the 
available data does not permit a reliable estimation of the 
power spectrum. In these cases, it is still possible to quan- 
tify the variability directly from the observed light curves by 
calculating the fractional rms variability amplitude (or sim- 
ply fractional rms), er rms , i.e., the square root of the variance 
of the light curve normalized by the mean flux over the pe- 
riod of ob servation after sub tracting the experimental noise 
(see, e.g., Nandraet al. 1997). Obtaining variances from X- 
ray light curves is a much less intensive observational task 
than obtaining reliable power spectra. As a result, the frac- 
tional rms variability has been measured for a relatively large 
sample of AGN for a variety of sampling patterns with ASCA, 
RXTE, and, more recently, with Chandra and XMM-Newton. 

Possible correlations between the square of the fractional 



rms and the X-ray luminosity have been studied in the past 
JNandra et alll997tlTurner et al.ll999HMarkowitz & E delson 
2001). More recently, the ava ilability of mass est imates for 
a large number of AGN (see IWoo & Urrvl 120021 for a re- 
cent compilation of ~ 300 masses) has allowed the study 
of correlations between a^ a and bl ack-hole mass, Mrh 
JLu & Yul l200!t iBian & Zhad l200l iMarkowitz & E delson 
12004110'Neill et alJ2005l) . The fact that larger black holes are 
less variable for a given observation period suggests that the 
observed (anti)correlation is partially due to a corresponding 
difference in the size of the X-ray emitting region. 

The fact that the fractional rms can be easily obtained 
from the observed X-ray light curves, and that the num- 
ber of AGN with measured variances is rapidly increas- 
ing, has motivated the theoretical modeling of the af ms - 
Mbh correlation. A fruitful approach to addressing this 
problem is to assume a functional form for the underlying 
power spectrum that depends on mass, P{y\ Mbh), and cal- 
culate the square of the fractional rms that one would ob- 
tain when observing according to some sampling pattern, 
i.e., cr 2 ms (Ai , At s , At^, Mbh)- Here, the set of values 
{At a , At s , Aib} stands for the duration of the observation, 
the sampling interval, and the binning interval respectively. 

The connection between the rms variability and the un- 
derlying power spectrum is crucial in theoretical models of 
the of ms - Mgjj correlation. As a first a pproxima tion, previ- 
ous works (Papadakis 2004; Nikolaiuk, Papadakis, & Czernv 
120041 but see also O'Neill et al. 2005) have assumed that the 
relationship between the observed fractional rms and the un- 
derlying (model) power spectrum is given by 

rl/2At s 

a r 2 ms (Ai , At s , Ai b , M BH ) = 2 / P{v; M BH ) dv , 

(1) 

where sampling effects are considered on the right hand side 
only through the finite limits of integration. To our knowl- 
edge, however, a systematic study addressing how good this 
approximation is for the typical sampling patterns employed 
in AGN observations has not been carried out yet. 

In this paper, motivated by current detailed timing studies 
of X-ray binaries, we consider model AGN band-limited noise 
spectra that can be described as a sum of Lorentzian compo- 
nents. This framework exposes (and allows to explore further) 
the underlying assumption often made when measuring AGN 
masses by comparing their power spectra with those of X-ray 
binaries, i.e., that the global timing properties of both systems 
simply "scale" with mass. An appealing feature of this frame- 
work is that it allows an analytical approach to account for 
the effects that the sampling pattern imprints on the observed 
power spectrum, facilitating the inclusion of sampling effects 
in theoretical models of the (7 2 ms - Mbh correlation. 

The paper is organized as follows. In Sj2]we discuss, from 
first principles, how the underlying power spectrum, whatever 
its functional form may be, differs from the observed power 
spectrum due to the finite and discrete nature of the observa- 
tions. In i|3]we motivate the use of Lorentzians to model the 
broad-band noise spectra of AGN and find an analytical ex- 
pression to account for sampling effects in this type of mod- 
els. In fj^we demonstrate the excellent agreement between 
our analytical results and Monte Carlo simulations designed 
to account for the distortions suffered by the underlying power 
spectral density due to sampling effects. Finally, in f|5|we dis- 
cuss the potential implications of neglecting sampling effects 
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in theoretical models when deriving masses from the af ms - 
Mbh correlation. For convenience, some of the mathematical 
details are presented in the appendices. 

While writing this paper, we became aware of a calcu- 
lat ion similar to the on e that we present in carried out 
bv lO'Neill et alJ J2005f) . Their derivation accounts for sam- 
pling effects directly in the power spectrum as they would be 
present due to discrete sampling in a "continuous" monitoring 
campaign (i.e., At s = Atb). The derivation that we present 
here is more general and shows explicitly how the different 
processes which one subjects the light curve to (i.e., binning, 
sampling, and segmentation), affect the observed power spec- 
trum for a sampling pattern characterized by the set of values 
{At , At s , Athl- In this sense , our derivation is more along 
the lines of Ivan der Klisl 111 989). explicitly showing what hap- 
pens to the real light curve in the process of obtaining the data. 

2. SAMPLING EFFECTS IN THE POWER SPECTRUM 

From the theoretical point of view, it is common to model 
and characterize the spectral content of a variable stochastic 
process by specifying the underlying (or model) power spec- 
trum, P(v), as a continuous function of the frequency v. This 
spectrum can be thought of as associated with a process, C(t), 
that varies over time in a continuous way. In the context of 
X-ray variability, C(t) could be given by the "underlying" 
(as opposed to observed) continuous and infinite X-ray light 
curve. In principle, complete knowledge of the observable 
C(t) implies complete knowledge of the underlying spec- 
trum, P(v), and vice versa. Therefore, if we had unlimited 
access to the X-ray light curve, we could reconstruct the un- 
derlying power spectrum accurately for every frequency. 

In practice, however, we do not have access to the contin- 
uous light curve, C[t), but rather to a finite and discrete set 
of values {C (t„)} that results from observing the source for 
a period of time At Q every t n = nAt s , where At s is the 
inverse of the sampling rate and n = 0, N — 1 (without 
loss of generality we will assume here that N = At Q /At s 
is even). In order to work with reasonably good signal-to- 
noise ratios, the data are usually binned for an interval of time 
Atb < At s (objects are said to be monitored "continuously" 
when At s = Atb)- In what follows we will refer to the set of 
values {At Q , At s , Atb} as the "sampling pattern". 

The finite and discrete nature of the observations not only 
dictates that we only have access to a finite and discrete set 
of values of the observed power spectrum, i.e., {Po^k)}, 
with v k = k/At a and k = -N/2 - 1, ...,N/2, but also that 
through the effects of red noise leakage and aliasing ( D eeming ! 

van der Klis 1989), these values do not coincide with 
the values of the underlying (model) power spectrum at those 
frequencies (i.e., {P Q (vk)} ^ {P{ v k)})- Moreover, because 
of the stochastic nature of the underlying process, the ob- 
served values {P (vk)} differ, in general, from observation 
to observation. For a stationary process, however, the aver- 
age taken over different observations (with the same sampling 
pattern), is well defined and it can be used as a consistent 
estimator of the underlying power spectrum iPriestlevlll989t 
iPress et afll 19921) . The aim of this section is to find an equa- 
tion that relates the continuous underlying power spectrum, 
P{v), with the set of values {(P Q {vk))} obtained by averag- 
ing several observations with the same sampling pattern. 

Because our goal is to find an equation that involves the 
model power spectrum, P(v), we need to work with contin- 
uous functions. In order to do so, we describe the set of ob- 



served values {C (t n )} in terms of the continuous variable t, 
representing the time, as follows 

C (t) = C bm {t;At h )W fin (t;At ,At s ) . (2) 

In this expression, 

C bm {t; At h ) = — / C(t') W(t - t'; At h ) dt 1 (3) 

is the "continuous binned light curve", 

oo 

W fin (t; At , At s ) = W(t; At ) At s J2 S{t - nAt s ) (4) 

n— — oo 

is the "finite comb function", and W(t, At) is the boxcar 
function of width At given by 

tit/, a j-\ / 1 if- At/2 < t < At/2 
W(t,At) = { Q otherw / se " " 

We can best interpret the right hand side of equation (0 
piecewise. The convolution of the continuous light curve 
C(t') with a boxcar function of width Atb retrieves a new 
continuous function of time, Cbi n (t; Atb). At any given time 
t, the value Cbi n (t; Atb) is the average number of counts re- 
ceived in the time interval |t — f'| < Atb. The finite se- 
ries of equidistant delta functions ensures that the observed 
light curve is not zero only when we are sampling (i.e., when 
t is some multiple of At s inside the range of observation 
[— At /2, Ato/2]) 1 . Equation can be seen as the mathe- 
matical representation of the process to which we subject the 
real (continuous) light curve when we observe it. Defined in 
this way, C (t) is a function of a continuous variable, it is 
zero for all t ^ t„, and its integral over time is equal to the 
total number of counts detected during the the entire period of 
observation At D . 

In order to understand how the observed values {{Poi^k))} 
differ from the underlying power spectrum due to sampling 
effects, we need to calculate the Fourier Transform of the light 
curve affected by the finite and discrete nature of the observa- 
tions, C (t), given by equation 0. By virtue of the convolu- 
tion theorem, its Fourier Transform, C (v), is given by 

/oo 
CunW- At b ) W fin (u - v 1 - At G , At s ) dv' , (5) 
-oo 

where Gun(y\ Atb) and Wfi n (v; At Q , At s ) are the Fourier 
Transforms of the continuous functions Cbi n (t; Atb) and 
Wfi n (t; At Q , At s ), respectively. The Fourier Transform of 
Cbin(t; Atb) can be computed using the inverse of the convo- 
lution theorem as 

C bin (v; At b ) = -^-W(v; At b ) C[y) , (6) 
At b 

where C{v) is the Fourier Transform of the continuous light 
curve, C(t), and 

W(v- At b ) = Sm(7wAtb) = At b sinc(™At b ) (7) 

TTV 

is the Fourier Transform of the boxcar function of width Atb- 

In order to have statistically independent measurements, we must have 
At s > Att,, while, in order to minimize aliasing effects, we would like 
to have Ai s = At b . In practice, however, with the exception of short- 
term observations, the sampling time employed in AGN observations is much 
larger than the binning time, i.e., At B S> Atj,. 
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Using again the convolution theorem, we can write 
the Fourier Transform of the finite comb function, 
Wfi n (t; At , At s ), defined in equation (0} as 

/oo 
E ff (i/; At s ) W{v - v'- Ai ) dv' . 
-OO 

(8) 

Here, W(v; At ) is the Fourier Transform of the boxcar func- 
tion of width At , W(t; At a ) (see eq. Q), and S^O; At s ) 
stands for the Fourier Transform of (At s times) the "infinite 
comb function" at equidistant times (see equation 0) 



oo 



E 

m— — oo 



S 5 (^;Ai s ) = > ' 5\ u- 



in 

At, 



(9) 



which is an "i nfinite comb function" at equidistant frequen- 
cies (see, e.g., Morse & Feshbach 1953). This expression al- 
lows us to carry out the integral in equation (|8} which yields 

oc 

Wf in {v\ At Q , At s ) = At D sinc[n(u - v m )At ] , (10) 



ra— — oo 



where we have defined the frequencies u m = m/ At s . 

Finally, using equations l|6}, @, and d 1 Oft we can write the 
Fourier Transform of the function C (t) in equation (0 as 



C {v) = £ C*wa(^ 



(11) 



where 



Cwa(^ — v m ) — At a / smc(nv' At h ) C(v') 

J — oo 

x sinc[7r(i/ — v 1 — v m )At ] dv' , (12) 

can be interpreted as the "windowed" and "aliased" contin- 
uous Fourier Transform associated with the real light curve 
C(t). Because of sampling effects, C (v) is a periodic 
function of period 1/At s (i.e., twice the Nyquist frequency, 
i^Nyq = 1/2 At s ) and thus the frequency v must be chosen to 
vary over a frequency range spanning at most 1 / At s . It is cus- 
tomary to choose this range as [— fN yq , ^NyJ- Even when this 
is the case, however, frequencies outside this range will influ- 
ence the value C (v) inside [— i^sr yq , ^Nyq]- When evaluated 
at the Fourier frequencies Vk = k/At Q , the set {C (vk)} is re- 
lated to the observed set of discrete Fourier coefficients {c&} 
that we would have obtained by taking the discrete Fourier 
Transform of the set of values {C (t n )} via 



C (v k ) = At a c k = £ C Q {t n ) e 



N-l 



N 



(13) 



The set of equations Q Q and (I12> show explicitly how the 
discrete set of values {C (vk)} is related to the underlying 
spectral content of the variable process, C{v). 

The relationship between the power spectrum affected by 
sampling effects and the Fourier Transform of any given ob- 
served light curve is given by P Q {v) = \C (v)\ 2 / At 2 . There- 
fore, using equation ll 1> . we can write 



Ai 2 £ ^WA(^ - Vm) C'wA^ ~ v m') >(14) 



where C^pJ^v) stands for the complex conjugate of Cwa^)- 
However, when the underlying process is stochastic, the vari- 
ance associated with t his estimator obtained from any given 
observation is 100% (lPriestlevlll989l iPress et alJl 19921) . In 
practice, a consistent estimator of the underlying power spec- 
trum is obtained by considering the average value of the ob- 
served power spectra obtained from many observations with 
the same sampling pattern, {At , At s , Ath}- For a stochastic 
process, the average over observation s of the cross-ter ms, i.e., 
m' 7^ m, in equation (I14t vanishes (Priestlev 1989), so we 
can write it as 



1 oo 



(15) 



1TL— — OQ 



Using equation (I12i . the average over observations becomes 

<^°M>=£ / / <CV)C>")) (16) 

x sinc(7iV A^b) sinc[7r(i^ — v' — v m ) At Q ] dv' 
x sinc(7ri/"Atb)sinc[7r(i' — v" — v m )At ] dv" . 

Furthermore, if the stochastic stationary process is ergodic, 
we can evaluate the average, (C(v') C*(v")), as an average 
over the ensemble of realizations {C(v)}. In AppendixlAlwe 
provide a brief demonstration showing that, in this case, 



<C(i/) C*{v")) = P(v') 6(v' - v") . 



(17) 



In this way, we can carry out one of the integrals in equation 
d!6t to obtain the relationship between the underlying power 
spectrum, P(v), and the estimator obtained from averaging 
over observations, (P (v)), as 

/>oo 

(Po(i>)> = / p W) W (v, v'; At , At s , At b ) dv' , (18) 



(19) 



where the function 

W (y, v; At , At h , At s ) = sinc 2 (7iVAi b ) 



x £^ sinc 2 [7r(^ — v' — v m )At ] 



ra— — oo 



is the "window function" that accounts for the effects of bin- 
ning, sampling and finite observation time, all at once. Strictly 
speaking, equation dl9l cannot be considered as the convolu- 
tion of the underlying power spectrum, P{v), with the func- 
tion VV [y, v'; At Q , At s , Atb) since the latter cannot be writ- 
ten as a function of the difference v — v' alone. 

As an aside, note that using the fact that the series of "sine" 
functions (weakly) converges to the Dirac delta, i.e., 



lim 77 sine (jxt)v) = 8{v) , 

rj — >oc 



(20) 



it is not hard to show that, as expected, in the limit of an in- 
finite and continuous observation the average of the observed 
power spectra tends to the underlying power spectrum, i.e., 



lim ( lim 

Ai Q — >oo \At B ,At b 



Ato<P (f)) = Piy) 



(21) 



m.m — — 00 



Whenever a sampling pattern is applied, however, power at 
any given frequency will "leak" to other frequencies, and vice 
versa. The extent of this spectral leakage is fully characterized 
by the function VV {v, v'; At Q , At s , At h ) (see Fig.[T}- 



Sampling Effects in the of ms - Mbh Correlation 



5 



» 




In many cases, it has been possible to obtain very good 
fits to the observed power spectra (of both neutron stars 
and black holes) using a small set of Lorentzians. It 
is important to point out that, whenever checked against 
each other, Lorentzian models usually yield better fits than 
power laws. Even when the fits are similarly good, the 
use of Lorentzians allows us to "follow" over time the 
different characteristic frequencies appearing in the power 
spectrum. This has been key to uncovering several cor- 
relations among the differe nt broad-band noise c ompo- 
nents in X-ray binaries ( Psaltis, Belloni, & van derKlisll 999: 
Belloni, Psaltis . & van der Klisl2002l) . 

A typical Lorentzian model for the underlying power spec- 
trum of an X-ray binary can be written as 



N L 

1=1 



Nl 2 

i=i 



A, 



A? + (i/-*f) 



■0\2 



(23) 



where Nl is the number of components, and rj, A/, and 
stand for the fractional rms, the width, and the centroid fre- 
quency corresponding to each Lorentzian. Defined in this 
way, the power spectrum is "two sided" and therefore 



F IG. 1 . — Grey-scale representation of the "window function", W (see 
eq. 1191 1. as a function of v and v' for At s = Atj, = 1024s and At = 
1024 2 s. Lighter shades of grey indicate higher values of VV. The function VV 
reaches its maximum at u' = with subsidiary peaks of decaying strength 
on both sides. As expected, the function VV is periodic in v with period Ai s . 



3. LORENTZIAN MODELS FOR AGN VARIABILITY 

In the previous section we showed how the average power 
spectrum inferred from a set of observations with the same 
sampling pattern differs from the underlying power spectrum. 
We can now write the average fractional rms that we would 
measure from such a set of observations as 



Li(v) dv . 



(24) 



(<7 r 2 ms (Ai , At Sl Ai b ; Mbh)) = 

N/2-1 



Although the data currently available for AGN might not 
allow yet for a clear distinction between power-law models 
and Lorentzian models, the underlying assumption of simi- 
lar X-ray variability properties between AGN and X-ray bina- 
ries advocates for the development of a common phenomeno- 
logical framework to study both kind of sources on the same 
footing. Motivated by this, we propose to use a small set of 
Lorentzians to model AGN power spectra. 

3.2. Sampling Effects in Lorentzian Models 

Broad-band noise components are usually characterized by 
ratios A;/z/ z ° 3> 1 and, in most cases, it is a good approx- 



2 y } j (P (z/fc; Mbh)) + (-PoC^Nyqi -Mbit)) ,(22) imation to set = in equation d23l and describe the 



k=l 

where the function (P (z/)) is the power spectrum affected by 
sampling effects and is related to the underlyi ng p ower spec- 
trum, P(y), according t o eq uatio ns (|18> and (1191 . Equation 
\22\ (together with eqs. 1181 and 1191 ) must be at the core of 
any theoretical model addressing the functional form of the 
°rUs~ Mbh correlation and considering the distorting effects 
imprinted in the data by the particular samplin g pattern used 
to obtain the observed fractional rms (see also. lO'Neill et al.l 
2005, where the case Aib = Ai s is discussed). 



3.1. Advantages of Lorentzian Models 

Inspired by the similarities shared by the the global timing 
properties of AGN and early observations of X-ray binaries, 
previous studies involving models of broad-band noise AGN 
power spectra are mostly based on underlying power spectra 
characterized by broken power laws. However, over the last 
decade, thanks to the timing capabilities of RXTE, it has been 
possible to study in great detail the X-ray variability of many 
galactic sources. It is now clear that the power spectra of X- 
ray binaries exhibit a rich morphology which, in general, can- 
not be adequately described in terms of broken power laws. 



band-limited n oise as a sum of zero -centered Lorentzians 
(Belloni. Psaltis. & van der Klis 2002). Let us then assume 
an underlying broad-band noise power spectrum that can be 
described by a set of Nl broad Lorentzians 2 , 



1 



7T * — ' 

{=1 



V? + V 2 



(25) 



where vi is the HWHM of the I th zero-centered Lorentzian 
but also plays the role of a "peak frequency" in the vP{y) vs. 
v representation (see Fig. [2}- 

With the normalization chosen in equation ( 1251 . the vari- 
ance of the light curves associated with the underlying power 
spectrum, P{v), is indeed the fractional rms. In other words, 
the total fractional rms corresponding to the band-limited 
noise spectrum is just 



2 

rms 



5 N L 

P{v)dv = Y,rt 
° 1=1 



(26) 



2 If required, the present formalism, including the derivation in Appendix 
IbI can be generalized straightforwardly to include narrow Lorentzian fea- 
tures, like the ones used to describe quasi-periodic oscillations. 
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TABLE 1 

Sample Lorentzian Model for the Power Spectrum 



NOTE. — Parameters defining the sample power spec- 
trum as a sum of three zero-centered broad Lorentzians. The 
fractional rms and the ratios between centroid frequencies for 
each Lo rentzian are similar to those observed in Cyg X-l 
iPottschmidt et al. 2003), while the peak frequency of the low- 
est Lorentzian (i.e., vi ) is similar to the lo w frequency break 
observed in NGC 3783 i Markowitz et alJ2()()3l) . 



This choice, together with the fact that the fractional rms is di- 
mensionless, determines the units of all the quantities related 
to the power spectrum. In particular, the underlying (continu- 
ous) power spectrum has units of inverse frequency while its 
discrete counterpart (affected by sampling effects) is dimen- 
sionless. Because of this, the discrete values P {vk) must be 
multiplied by the observing time At Q in order to be compared 
with the underlying values P{yk)- This is also the reason for 
which the factor At Q appears explicitly in eq. (1211 . 

Figure |2 shows a sample Lorentzian model for the broad- 
band noise power spectrum of a hypothetical AGN (see also 
Tab.^. In this case, the fractional rms and the ratios between 
centroid frequencies for each Lorentzian are similar to those 
observed in Cyg X-l (Pottsch midt et al.l2003i) . while the peak 
frequency of the lowest Lorentzian (i.e., ) is similar to t he 
low frequency break of NGC 3783 (Markow itz et al.l2003l) . 

Using equations dl 8I > and dl9> and the model power spec- 
trum, P{y), given by equation (|25> we can write the average 
observed power spectrum, (P (z/)), as 



N L oo 



(27) 



I— I m— — oo 



with 



sine 2 fa„, (V) — x'] 



where we have defined the dimensionless variables 



dx' 



(28) 



x = TtAt v . 



A = At h /At , 

a m {v) = TiAt (v - v m ) 



This integral can be computed analytically in the complex 
plane using a slight modification of the residues theorem (see 
Appendix[E]for the details). The result is given by 



h,m{v) - 



4m( 
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A 2 2^ A i 

1=1 



Mi 

1 - (1 - A) cos(2Aa m ) 



(29) 
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FIG. 2. — The underlying power spectrum for a hypothetical AGN (solid 
line) modeled as a sum of three broad Lorentzians (dotted lines). The arrows 
indicate the positions of the peak frequencies (see Tab. 0. Vertical lines 
indicate the frequencies corresponding to the different timescales involved in 
the numerical experiment performed in J4] 



a m sin(2a ?i 



.-26, 



A 5 f..\_ (h ~ a m) „-Abi 



1 sinh 2 (A& ; ) 



2(al + b?)H? 
A 



sinh(A6 ; ) , 



2n 2 b 2 ' 



where we have omitted the frequency dependence of the co- 
efficients a m (v) on the right hand sides of the expressions for 
A\ m (y). The symmetry property 



A. 



I,-m(-") =A,rM 



(30) 



reflects the fact that sampling effects do not affect the parity 
of the observed power spectrum, i.e., (P Q (v)) = (P (—v)). 

It is important to summarize here what we have accom- 
plished with equations \21\ and J29l l. We have reduced the 
problem of accounting for the sampling effects caused by an 
arbitrary sampling pattern, {At , At s , Atb}, in any broad- 
band noise power spectrum that can be decomposed into a 
sum of zero-centered Lorentzians, to evaluating a simple sum. 
In practice, the sum in equation d27t converges extremely fast 
and it is necessary to add only a few (order of five) terms at 
any given frequency. This makes the process of accounting 
for sampling effects in Lorentzian models highly efficient if 
the aim is to run a grid of models to contrast them against 
observations in order to constrain the functional form of the 
underlying power spectrum. We stress that, if the underly- 
ing power spectrum is given by equation d25), the average 
(P (v)} calculated according to equations J27l > and d29t pro- 
vides the value to which Monte Carlo simulations, devised to 
account for the sampling effects imprinted by the sampling 
pattern {At Q , At s , Att>}, will converge. In the next section 
we show that this is indeed the case. 

4. CONTRASTING THE ANALYTICAL PREDICTION TO 
MONTE CARLO SIMULATIONS 

The aim of this section is two-fold. First, we want to test 
how well our analytical prediction is able to account for the 
distortions suffered by the underlying power spectrum due to 
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sampling effects. Second, we would like to obtain an idea 
about how important of an effect one is neglecting when cal- 
culating the rms variability using equation ([0 instead of i22\ 
for a sampling pattern typical of AGN observations. In order 
to do this, we performed a series of numerical experiments 
calculating independently both sides of equation i22\ assum- 
ing the model power spectrum characterized by the three 
broad Lorentzians shown in Table^see also Fig.|2j. 

In order to test the validity of equation i22\ . we need 
to provide a good representation of the 'infinite" and "con- 
tinuous" light curve, C(t), associated with the underlying 
power spectrum, P{v). Then we can segment, bin, and sam- 
ple this "real" light curve according to the sampling pattern 
{At , At s , A^b}, obtain a set of time series {C (t n )}, and 
calculate the rms amplitude, rr 2 ms , associated with each sam- 
pled light curve as 



N-l 



N- 



Co 



(31) 



where N = At / At s and Co is the average corresponding to 
the particular time series {C (t n )}. The mean value in equa- 
tion J22L (<7 2 ms ), can then obtained by averaging over several 
"observed" light curves with the same sampling pattern. 

Of course, a numerical representation of the "real" light 
curve C(t) will still be given by a discrete set of values 
{C(ti)}, much in the same way as we can only evaluate the 
underlying (continuous) power spectrum at a finite number of 
frequencies {P(yk)}- The way of doing this is to consider a 
minimum frequency, v m i n , and a maximum frequency, v max , 
such that the contributions to the rms variability from power 
outside the range [v m i n , u max ] is negligible, i.e., 



P{u)dv, I P(v)dv < 



P[y) dv . (32) 



In the case of the model under consideration this can be ac- 
complished by choosing v mm <C v\ and i/ max 3> v%. Note 
that this is a much stronger condition than requiring that 
fmin -C 1/At and/or i/ max 3> 1/2 At s for the particular 
sampling pattern under consideration. These conditions over 
Vmin and v max are necessary but not sufficient. 

For all practical purposes, the time series {C(ti)} obtained 
in this way can be considered as an accurate representation of 
C(t) in the sense that the power spectrum estimated from it 
will not be noticeably affected by sampling effects. Note that, 
the set of values {P (yi~)} corresponds to one given realiza- 
tion of the stochastic process and, therefore, in general, it does 
not coincide with the values of the underlying power spectrum 
{P{vkY\ at any given frequency even in this case! However, 
the lack of sampling effects is noticed in the absence of the 
trends that would be introduced by red noise leakage and/or 
aliasing inside the range [z/ min , u max ]. 

In order to set up the numerical experiment, we generate 
the "continuous" light curve by setting ^ m j n = 1/ AT max = 
9.3 x lCr 10 Hz (i.e., AT max = 2 30 s or roughly 36 yrs) and 
*W = ^Ny q = l/(2AT min ) = 7.8 x l(T 3 Hz (AT min = 
64s). This yields N = AT max /AT min = 2 24 ~ 10 7 ini- 
tial "data" points. The frequencies corresponding to v m \ n and 
'-'max are indicated with vertical lines (dotted and triple-dot- 
dashed line, respectively) in Figure |2] The conditions stated 
in equation d32i are clearly satisfied. 

Under our current set of assumptions, we can generate the 
time series {C(t;)}, with U = lAT min and I = 0, N - 1 
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power spectrum, P(v), as a function of the observation time At calculated 
in three different ways. Filled circles: results obtained from Monte Carlo 
simul ation s. Solid line: the analytical prediction given by the right hand side 
of eq. 1221. D ashed line: results obtained by neglecting sampling effects, i.e., 
using ea. 111. In this case, the binning time is equal to the sampling time. 



from the set {P(v k )}, with v k = k/AT max and k = -N/2 + 
l,...,N/2, as 



N/2 

C(ti) = E Ck ' 

fc=JV/2-l 



with 



Ck = 



P{y k ) 



D 27riip(u k ) 



AT„ 



(33) 



(34) 



In order to generate a true realization of the stochastic pro- 
cess with underlying power spectrum P{v), the amplitudes 
under the square root in equation (I34> must satisfy P(vk) — 
P(—Vk) and be drawn from a x 2 distribution with two de- 
grees of freedom for all ^ ^Nyq and from a \ 2 distri- 
bution with one degree of freedom for i^Nyq- Alternatively, 
this can also be done by generating two normally distributed 
random variables, multiplying them by [P(ffc)/AT m ax] 1 / 2 , 
and using them to define the real and imaginary parts of the 
discrete Fourier Transf orm given by the set {cfe}(see, e.g., 
Timmer & Konig 1995, for more details). Moreover, in or- 
der to ensure that the generated light curve is real, the set of 
random phases {<j){vk)} must satisfy {4>{— Vk)} = {~0(^fe)} 
for all k + 0, Nj2 and </»(0) = or tt, as well as <X^N yq ) = 

or TT. 

Note that, by construction, if we were to take the dis- 
crete Fourier Transform of the set {C(ti)} and calculate 
{P(vk)} = AT max {|cfe| 2 } then sampling effects would be 
negligible in the range [t-Wn, v max ]. For all practical purposes 
then, the set {C(t;)} can be thought of as an accurate repre- 
sentation of the "infinite" and "continuous" light curve C(t). 

We can now define a sampling pattern, {At D , At s , Aib}, 
and do everything we would do if we were recording a time 
series from the continuous light curve C[t). In order to be 
able to compute an average value for the rms variability we 
divided the original light curve of duration AT max in 16 light 
curves of duration Ai . max = 2 26 s. We also defined the 
duration of each time bin as A<b = 1024s, as a representa- 
tive value used in AGN observations. Note that this value of 



8 



PESSAH 



TABLE 2 

Sampling Effects on the Fractional RMS 



n 


At /At s 


2 / 2 
<J rms,WA/ <J rms,NS 


3 


8192 


~ 10.00 % 


4 


4096 


~ 20.00 % 


9 


128.0 


~ 100.0% 


12 


16.00 


~ 200.0% 



NOTE. — Comparison between the fractional rms ob- 
tained considering and neglecting sampling effects as a 
function of the ratio between the observing and the sam- 
pling times. The symbols of mS]WA and of ms NS de- 
note the fractional rms derived f rom the "windowed" and 
"aliased" power spectrum (eq. 1221 1. and the value ob- 
tained when sampling effects are neglected (eq. HI). 



Atb will include 16 discrete points of the original time series 
{C(ti}}. To simplify the discussion, in what follows, we will 
consider a continuous monitoring campaign, i.e., At s = Atb 
(we consider a case in which At s > Atb in the next section). 

We calculated the fractional rms variability for each light 
curve of duration At Q _ max according to equation (13 1 1 and ob- 
tained the average value that we took as the left hand side of 
equation i22\ . We repeated this procedure defining an increas- 
ingly larger number of segmented (independent) light curves 



of duration At = 2~™Ai 



with n = 0, ...,13. Fig- 



ure|3]shows the values of (of ms ) as a function of the observing 
time 3 , At . The analytical prediction (i.e., the right hand side 
of eq. 1221 ) is shown with a solid line, the filled circles repre- 
sent the average values (Crms) obtained from the light curves 
generated with the Monte Carlo simulations, and the dashed 
line shows the value of of ms as obtained using equation 

It is important to stress that we have not adjusted the nor- 
malization of any quantity when plotting the different vari- 
ances in Figurey] The convergence of the different fractional 
amplitudes for large values of At to the total rms amplitude 
of the model (i.e., of ms = 0.09) indicates that sampling ef- 
fects become negligible when the harmonic content of the ob- 
served light curve is a good representation of the underlying 
power spectrum. This is in agreement with equation (12 1 i . As 
the observation time decreases, the fractional rms, as calcu- 
lated according to equation 0, underestimates the values ob- 
tained from either side of equation i22\ . i.e., the fractional rms 
obtained from either the analytical result or the Monte Carlo 
simulations, by a factor of up to 200% (see Tab.|2]for details). 

In the previous numerical experiments, the fractional rms 
was obtained directly from the light curves without the need 
of deriving the power spectra associated with them. Fig- 
ure|4]shows the average power spectra derived from the Monte 
Carlo simulations (filled circles) for the case At s = Atb = 
1024s and At = 65536s (corresponding to n = 10) together 
with the corresponding analytical result (open circles) derived 
using equations (I27M291 . The distorting effects suffered by 
the power spectrum derived from the finite and discrete light 
curves with respect to the underlying power spectrum (solid 
line) are evident. This remarkable agreement between the av- 
eraged observed power spectra derived analytically and using 
the simulated light curves is the reason behind the coincidence 
of the associated fractional variability amplitudes in Figure[3] 

3 We have used the dimensionless quantity v\ At as the independent vari- 
able. In this case, sampling effects are not important for v\ At Q < 1. 
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FIG. 4. — Sampling distortions suffered by the model power spectrum 
defined in Tab.fT1t.so/it/ line), as obtained from calculating the (average) dis- 
crete power spectrum associated with the light curves generated via Monte 
Carlo simulations (filled circles), and as predicted by our analytical calcula- 
tion (open circles). The sampling pattern is defined by At = 65536s and 
At B = At]-, = 1024s. Note that the discrete values P (^fc) must be multi- 
plied by At in order to be compared with the underlying values P(u^). 



As a corollary, we note that, by adopting a Lorentzian 
mode l for the underlying power spectrum, the results from 
£13-21 can be used to test Monte Carlo simulations when the 
frequencies are assumed to be uncorrelated. 

5. DISCUSSION 

This is the first paper in a series aimed to exploit the sim- 
ilarities observed in the X-ray timing properties of AGN and 
X-ray binaries in order to measure the masses of supermassive 
black holes using X-ray variability. 

The idea of using X-ray variability to obtain the masses 
of AGN by comparing their timi ng properties to those of 
X-ray binaries is not new ffiavashida et al. 199 81 ILu & Yul 



2001 



Czernv et all 1200 It iBi an & Zh aol 12003b tPapa dakis 
2004; Nikolaiuk, Papada kis. & Czernvl2004l) . A key assump- 
tion necessary for this procedure to provide reliable results 
is that the processes driving the variability, and therefore the 
associated underlying power spectra, are similar. Although 
the data cur rently ava ilable for some AGN is encoura ging in 
this respect dUtflev et all2002tlMarkowitz et all2003l) . so far 
there have been no attempts of describing the timing proper- 
ties of both kind of objects on the same footing. Motivated by 
the latest timing studies carried out on a number of galactic 
black-hole candidates, we have proposed to model the broad- 
band noise spectra of AGN in the same way that has been suc- 
cessful with many X-ray binaries, i.e., as a sum of Lorentzian 
components. 

A particular concern when connecting theoretical models 
of broad-band noise X-ray variability with observations is that 
sampling effects (especially red noise leakage) are much more 
important for AGN than for X-ray binaries. In the last few 
years it has become evident that in order to derive reliable 
power-spectra estimates from the observatio ns it is impe rative 
to account for these sampling effects JlJttlev et alJ2002l) . Less 
attention has been paid, however, to the effects that sampling 
has on the connection between the underlying power spectrum 
and the rms variability (but see O'Neil let all2005l) . 

In order to illustrate the importance of this latter issue, 
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TABLE 3 

Typical Sampling Patterns for AGN 



Term Type 


At a 


At 8 




2 / 2 
%s,WA'%s,NS 


Long 


4.27x266 d 


4.27 d 


1024.8 s 


~ 78% 


Medium 


3.20x151 h 


3.20 h 


1152.0 s 


~ 32% 


Short 


2000 x84 s 


2000 s 


2000.0 s 


~ 77% 



NOTE. — The sampling patterns for short-term observations are 
typical of ASCA, Chandra, and/or XMM-Newton. The values of 
{At , At s , Atb} for long-term observations are more representative of 
RXTE observations. The meaning of the symbols f^ms WA an( ^ CT rms NS 
is the same as in Tab.l2l 



let us assume a hypothetical Lorentzian model resembling 
the underlying power spectrum obtained for NGC 3783. 
Markowitz et al. (2003) found that a reasonable description 
of this power spectrum can be obtained by a double bro- 
ken power-law model with characteristic ("low" and "high") 
break frequencies given by v\ — 2.00 x 10 _7 Hz and Vh = 
3.98 x 10 _6 Hz, and an amplitude A ~ 0.01 (the correspond- 
ing power-law slopes at low and intermediate frequencies 
were fixed to and -1 respectively and the best fit at high fre- 
quencies is consistent with a power-law of -2). A Lorentzian 
model to describe this underlying power spectrum could be 
obtained using two Lorentzian components (dotted lines in 
Fig. |5} as follows 



P{y) 



V-2 



7T uf + V 1 



(35) 



where we have denoted the characteristic frequencies as v\ 



vi and V2 = Vh- We assume here that rf 



0.075 so 



that the peak amplitude of each Lorentzian in vP(y) space 
is roughly 0.01. As an aside, we note that the ratio of fre- 
quencies derived from observations is such that these two 
Lorentzians (with the same amplitude) lead naturally to a fiat 
slope in the vP(v) representation. This feature is observed 
in several X-ray binaries, i.e., in many cases the frequencies 
and amplitudes of the different broad Lorentzians are such 
that the overall shape of the (underlying) power spectrum re- 
sembl es a 1/V power law over a wid e range in frequencies 
JBelloni. Psaltis. & van der Klisl2002l) . 

Let us now suppose that we "observe" the model defined 
by equation d35l according to the different sampling patterns 
involved in the various types (long-, medium- and short-term) 
of monitoring campaigns designed to observe NGC 3783 (see 
Tab. |3- By calculating the rms variability according to equa- 
tion Q, (of ms NS , i.e., not considering sampling effects), and 
equation ( 1221 . (of ms WA , i.e., when the rms variability is cal- 
culated based on the "windowed" and "aliased" power spec- 
trum) we can study how important are the consequences of 
neglecting sampling effects. The last column in Table[3]shows 
the ratio of these two variances for the different types of sam- 
pling patterns. For the adopted power spectrum, neglecting 
sampling effects can lead us to underestimate the observed 
value of the rms variability by up to 80%. These differences 
can also be understood by looking at the sampling effects di- 
rectly in the power spectrum. Figure |5] shows the "observed" 
power spectra that we would obtain from monitoring a source 
with an underlying power spectrum given by equation J35I . 
according to the sampling patterns listed in Table|3] 

It is not straightforward to see how these differences in the 
variances will affect the mass estimates derived from a theo- 
retical model of the of ms - Mbh correlation. This will depend, 
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FIG. 5. — Illustration of the distorting effects suffered by the underlying 
power spectrum, P(y) (so lid l ines), modeled with two Lorentzian compo- 
nents (dotted lines, see eq. I35J) when it is "observed" according to different 
sampling patterns (see Tab.131. The different symbols, filled diamonds, open 
circles, and filled triangles, indicate the "windowed" and "aliased" power 
spectra that would be derived from long-, medium-, and short-term observa- 
tions respectively. 



of course, upon the functional form of the of ms - Mbh corre- 
lation itself. We note here that a difference in rms variability 
will likely be amplified when translated into a difference in 
mass. This is because the observed fractional rms is a rather 
fiat function of black-hole mass. Indeed, for a large number 
of AGN the fractional rms is known to vary between 30% and 
40% in long-term observations and between 2% and 30% in 
short-term observati ons, over al most four orders of magnitude 
in black-hole mass (Mar kowitz et all2003l) . 

This flattening of the observed cr r 2 ms - Mbh correlation for 
larger observation times suggest that, for the sake of deriving 
masses, short observations should be favored over long term 
monitoring campaigns. This is, of course, when sampling ef- 
fects (in particular red noise leakage) will affect the most the 
connection between the model power spectrum and the ob- 
served fractional rms variability. For this reason, it is vital to 
incorporate sampling effects in theoretical models that aim to 
estimate black-hole masses via the of ms - Mbh correlation if 
they rely on a parametrization of the underlying AGN power 
spectrum in terms of Mbh ■ 

It is evident from Figure[5]that the effects of spectral leak- 
age have to be considered when (continuous) model power 
spectra are to be compared against (average) power spectra 
obtained from observed time series. State of the art techniques 
to "convolve" the models with the observation sampling pat- 
terns currently involve sophisticated Monte Carlo methods. 
These simulations are crucial once the "best fit model" has 
been found in order to assess the confidence levels associated 
with the best fit parameters. The high numerical costs asso- 
ciated with this type of simulations, however, usually rules 
in favor of mathematical models that can be described with a 
low number of free parameters which might not allow enough 
freedom to properly describe the data. 

The proposed Lorentzian models for broad-band noise vari- 
ability offer a tremendous advantage in this sense. In the mea- 
sure that the available data allows it, the gain in computational 
speed obtained by using analytical expressions to account for 
sampling effects could be used to explore parameter spaces of 
higher dimensions. This will play in favor of more realistic 
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APPENDIX 

AVERAGE OVER REALIZATIONS 

In f|3]we stated that the average over ensembles of the quantity {C{v')C* (v")) was related to the underlying power spectrum 
of a stationary stochastic process via equation d!7l >. Here, we provide a brief demonstration of this statement. 

By definition, the average {C{v')C*(v")) can be written in terms of the corresponding average in the time domain as 

/oo />oo 
/ (C(t')C*(t")) e-' 2 ™'*' e* W '*"W , (Al) 
-oo J — oo 

where C(t) corresponds to the inverse Fourier Transform of C(y)- For a continuous process with zero mean, the quantity 
(C(t')C* (t")) is the correlation function , t"). Furthermore, if the process is stationary, the expectation value (C(t')C* (t")) 
can only be a function of the time difference t' — t" and, therefore, 

/>oo />OC 

,l\r<* t,,U\\ / / t/V 4-ll\ —i2-!tv't' A2iru"t" Jj JjH 



{C(v')C*{v"))= I / — t") e ~ 1 e * dt'dt" . (A2) 

J —oo J — OO 

Defining r = t f — t" and f = we can write 

pOO poo 

{C{v')C*{v")) = / / £(T)e- l2 ™' f e l2 ™"( f - r Wdf , (A3) 

— oo J —oo 

pOO 

£(T) e -i***"T dr / e -i^'-"")r df , (A4) 



= F[t](v")5{v' -v"), (A5) 

where the quantity F[£](v") stands for the Fourier Transform of the autocorrelation function, which is the underlying power 
spectrum, P{v"). In this way, we obtain the result quoted in equation d!7l . i.e., 

{C{v')C*{v")) = P(y") 8{v' - v") . (A6) 
CONTOUR INTEGRAL 

In Sj3 .21 we presented an analytic expression for the coefficients Ii : m{v) involved in the calculation of the average power 
spectrum affected by sampling effects according to the sampling pattern { At , At s , A£t>}- Here, we outline how to obtain them. 

The integral defining the coefficients Ii. m (v), i.e., 
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fh f°° sin 2 (Ax') sm 2 [a m {is) -a/] dx' 



k ' m{v) ~ n (Ax') 2 [a rn { V )-x>Y bf + x' 2 ' (B1) 

exists and is real for any non-zero value of 6; and any finite value of A, and a m (v). A convenient way to calculate its value 
is to allow the variable x' to be complex and compute Ii t m{v) as a contour integral in the complex plane by choosing suitable 
integration paths. The idea, as usual, is to find a contour of integration that contains the real axis and close it conveniently in 
such a way that the imaginary contribution vanishes when the real part of the contour extends from — oo to oo. In order to do so, 
it is necessary to rewrite the sines in terms of co mple x exponentials that decay in the upper or lower half of the complex plane 

aation dB 1 i as 



depending on their sign. We can rewrite equation (IB It as 



p=i ± 

where the functions forp — 1,2,3 are given by 

f±(z) = {e 2iAz -l)e ±i M z - a ™y , 
f±(z) = (e- 2iAz -i) e ±mz-a m )] ; 

f±(z) = 2(l-e ±2iA *). (B3) 

Note that because A = At^/ At < 1, the functions /+ and f~, for p = 1,2, 3, decay exponentially in the upper and lower 
halves of the complex plane respectively. A minor complication arises because the integrand presents (simple) poles on the real 



Sampling Effects in the of ms - M B h Correlation 1 1 



▲ 




T 

FIG. B6. — The integration paths in the upper (F_)_) and lower (r_) halves of the complex plane used to calculate the coefficients /; m (^) in equation <BH . 
The functions fp defined in equation IB3I decay exponentially over T± for lm(z) ^ 0, respectively, when the radii of the contours tend to infinity. 

axis, i.e., at z = and z m — a m {v), so some care is required when choosing the contours of integration in either half. Figure lB6l 
shows a possible way of choosing two closed paths. The contours T± are adequate to calcula te the integrals cont aining the 
functions f^, respectively. Using a slight modification of the residues theorem (see, e.g.. lMarsden & Hoffmanlfl 999ft we can 
calculate the integrals in equation (IB2i to obtain 

2 b 3 — > 

IiM = ^ ±2*ines(£,±bi) ± nines(f±,Q) ± TriTZes[f±, a m {v)] , (B4) 

p=i ± 

where lZes(f, z) stands for the residue of the function / evaluated at the (simple) pole z. Note that the signs in this equation 
properly take care of the negati ve or ientation of the contour T_. After some lengthy, but otherwise straightforward, algebra to 
regroup similar terms, equation JB4> yields the result quoted in equation \29\ . 
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